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ABSTRACT 

The evolution of the spin parameter of dark halos and the dependence on the 
halo merging history in a set of dissipationless cosmological ACDM simulations is 
investigated. Special focus is placed on the differences of the two commonly used 
versions of the spin parameter, namely A = JE 1 ^ 2 G~ 1 M~ 5 / 2 (Peebles 1980) and 
A' = J/[V2A/ vir i? vir Kir] (Bullock et al. 2001). Though the distribution of the spin 

transfer rate Tj[ := Ag^/Aj^ t (which is the ratio of spin parameters after and prior to 
merger) is similar to a high degree for both A and A', we find considerable differences 
in the time evolution: while A' is roughly independent of redshift, A turns out to 
increase significantly with decreasing redshift. This distinct behaviour arises from 
small differences in the spin transfer distribution of accreted material. The evolution 
of the spin parameter is strongly coupled with the virial ratio r\ := 2i?kin/|-E'pot| of 
dark halos. Major mergers disturb halos and increase both their virial ratio r\ and spin 
parameter \W for 1 — 2 Gyrs. At high redshifts z = 2 — 3 many halos are disturbed 
with an average virial ratio of rj ~ 1.3 which approaches unity until z = 0. We find 
that the redshift evolution of the spin parameters is dominated by the huge number 
of minor mergers rather than the rare major merger events. 

Key words: cosmology: dark matter halos, angular momentum, spin parameter, 
mergers - methods: numerical 



1 INTRODUCTION 

It is among the puzzling tasks within the theory of galaxy 
formation to understand the origin and history of angular 
momentum both of dark matter halos and embedded galac- 
tic disks. One reason for the high relevance of galactic spin 
is that it determines the most basic properties of a galactic 
disc, its surface density distribution and its spatial exten- 
sion through the equilibrium of gravitational and centrifu- 
gal forces, as discussed already by Hoyle as early as 1949 
(Hoyle 1949). A kind of basic picture of galactic evolution 
has emerged in the undertow of the seminal works of White 
& Rees (1978) and Fall & Efstathiou (1980), which connects 
the angular momentum of a dark halo to its mass growth 
history (see also Blumenthal et al. 1986, Mo, Mao & White 
1998, van den Bosch 2000, Cole et al. 2000, Somerville, Pri- 
mack & Faber 2001, Kauffmann, White & Guiderdoni 1993). 
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Until the beginning of this century it was considered as 
established that bound cosmological structures gain most 
of their angular momentum by the tidal effect of the large 
scale mass distribution: After a merger the remnant is left 
with the angular momentum imprint of the tidal field. In 
recent years it was realized that yet another effect is con- 
tributing to a halo's angular momentum history to at least 
the same degree as large scale tidal forces, namely its dis- 
crete mass acquisition history (Gardner 2001, Vitvitska et al. 
2001, Mailer, Dekel & Somerville 2002). In this scenario, or- 
bital angular momentum is transferred to the remnant's in- 
ternal angular momentum during a merger. Statistically, the 
net angular momentum of a large number of mergers would 
be zero (if infall occurs from random directions). However, 
due to the low number of major mergers that occur during a 
halo lifetime, randomization is ineffective and there remains 
an imprint on the halo's spin parameter by the final (and 
often only) major merger. Mergers of dark matter halos and 
galaxies thus moved to the center of interest not only in the 
field of early type galaxy formation but also in the context of 
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angular momentum acquisition (Burkert & D'Onghia 2004. 
D'Onghia & Burkert 2004). 

The angular momentum of dark matter halos is usually 
parametrised in terms of a dimensionless quantity, the spin 
parameter A. This parameter appears in two versions; the 
first, original one was introduced by Peebles (1980). A re- 
vised formula was later proposed by Bullock et al. (2001) 
in order to overcome the energy dependence of the classic 
version. 

In no study so far the differences between the classical 
A and the revised version A' have been investigated in detail. 
This is the main ambition of our work. For this purpose we 
analyse the data of 82 dissipationless numerical simulations 
of dark matter halo formation in a comoving cube of 20 
Mpc starting with redshift z=2.8. In order to understand 
our results we also investigate the connection between the 
spin parameters and the time dependent virial coefficient 
r) = 2E kiTi /\E pot \. 

In §2 we describe the numerical simulations in detail. §3 
analyses of the merger anatomy and statistics. A detailed 
investigation of the spin transfer in merger and accretion 
events is given in §4 and §5, with special attention drawn to 
the differences between the two A-versions. In §6 we investi- 
gate the virial parameter of dark halos and it's connection 
with the spin parameters. The global time evolution of A 
and A' is shown in§7. Finally we summarize and discuss our 
results in §8. 



2 NUMERICAL SIMULATIONS 

We employ the GRAFIC code which is part of the COSMICS 
package developed by E. Bertschinger (Ma & Bertschinger 
1995) to set up random Gaussian initial conditions. In to- 
tal, 82 dissipationless cosmological ACDM simulations arc 
performed with different random seeds for the initial fluc- 
tuation field, using a N-Body code with a Barnes-Hut tree 
method (Barnes & Hut 1986) in combination with the spe- 
cial purpose hardware GRAPES (Makino & Funato 1993, 
Ebisuzaki et al. 1993). Each run represents a 20Mpc co- 
moving cube with the standard cosmological parameters 
n m = 0.7, ft A = 0.3, h = 65 and cr 8 = 1.13 1 . The sim- 
ulations start at redshifts 40 < z ini t < 53 (depending on 
when the maximum overdensity reaches unity) and follow 
the nonlinear evolution of 64 3 particles until z = 0. Position 
and velocity data are stored at 21 different redshifts mar- 
king approximately equal time intervals. The gravitational 
potential is smoothed within 7.8 comoving kpc around indi- 
vidual particles, each of which represents 1.07 x 10 9 M Q of 
cold dark matter. 

We identify halos using the friends-of-friends (FOF, see 
e.g. Gotz, Huchra & Brandenberger 1998) algorithm with 
a linking length parameter of 0.2 in units of mean particle 
distance. Each FOF cluster of particles is truncated at a 

1 Note that these simulations were performed earlier than the 
current WMAP results (ct 8 = 0.9 ± 0.1, see Spergel et al. 2003) 
were published. Our value for erg is following the expression 
<t g = (0.5 . . . 0.6)Q-°- 56 from White, Efstathiou & Frcnk (1993). 
We don't expect this discrepancy to affect our results. Especially 
results concerning individual merger events should be completely 
unaffected. 



radius R v i T which defines a sphere around the most bound 
particle within which the mean density is A(z)p c (z), p c (z) 
being the critical density at redshift z. Following the spher- 
ical collapse model, A(z) would be constant in an Einstein- 
de Sitter universe, with A = 18-7T 2 ~ 178 (see e.g. Peebles 
1980). Due to the model dependence of the density param- 
eter Q(z) = (fi m — l) 3 r(z; Sl m , SIa) this does not hold for 
a low density model like the one considered here. Neglect- 
ing radiative contribution to the mean cosmological density, 
A(z) takes the general form 

A(z) = 18it 2 +82f(z)-39f(zf, (1) 
where 

_ n m {i + z f _ 
nz> n m (i + *)3 + n A [2) 

(Bryan k Norman 1998). 

Our simulations yield a total number of 10563 halos, 
each with more than 100 particles at z — 0. In order to avoid 
a mass truncation bias in our merger trees, we consider only 
halos which contain at least 500 particles at present epoch, 
corresponding to a minimum halo mass of 5.35 x 10 n M o 
which leads to a sample of 2383 halos; particle clusters as 
small as 50 particles (1.07 x 10 11 Mq) are nevertheless taken 
into account as progenitors and satellites at higher redshifts. 
Note that throughout this study we didn't detect any mass 
dependence of our results, which therefore should be valid 
for lower mass halos as well. 

We construct merger trees for the complete sample of 
halos with M(z = 0) > 5.35 x 10 n M Q in the following way: 
A major merger (MM) is defined to be an event where two 
initially separated halos appear FOF-bound in the next data 
output with 3 restrictions: 

• the less massiv, infalling clump (denoted as the satel- 
lite) has a fraction of at least van > 1/5 of the larger clump's 
(the halo's) mass. 

• the merger remnant contains more than 60% of the ini- 
tial satellite's particle reservoir 

• the massive progenitor halo has a mass of at least 40% 
of its mass at z — 0. 

The first restriction separates mergers from events which 
we will denote as "accretion" or "minor mergers" (mM) 
throughout. Note that we do not distinguish between mi- 
nor mergers and accretion; both terms correspond to any 
detected matter acquisition which increases the progenitor's 
mass by less than 20%. The second criterion guarantees that 
the mass growth after the preceding time step does in fact 
originate from the satellite infall and not from accretion or 
mM events which may well happen during the same time 
period. The third point finally rules out mergers which, due 
to the low mass of the remnants, do not have a significant 
impact on the dynamics of the final halo or galaxy at z — 0. 
It is clear that our particular choice of the parameter val- 
ues (60% and 40% in the 2nd and 3rd item, respectively) 
is somewhat arbitrary and lacks a quantitative justification. 
Nevertheless it seems plausible that this choice of numbers 
is suitable for the purposes described above. By identifying 
all progenitor-remnant relations throughout all data outputs 
we construct a list of progenitors of each halo at z = as 
well as the according mass acquisition history. 
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3 MERGER STATISTICS AND ORBITAL 
PARAMETER 

Before we are going to examine the spin evolution of dark 
matter halos in detail we want to "set the stage" by briefly 
discussing the abundance and the orbital characteristics of 
major mergers. 



3.1 Orbital Parameters 

Certainly the main quantity that characterizes the dynam- 
ics of a major merger is the mass ratio of the involved halos. 
The dynamical properties of the remnant, however, depend 
on both the internal and orbital energy and angular mo- 
mentum contributions. If we describe the approach of two 
halos following the Kepler laws for classical two body sys- 
tems (which is well justified as long as the objects are spa- 
tially separated) , the orbital energy and angular momentum 
content of the system is fully described by their orbital pa- 
rameters. The eccentricity is usually defined by 

1 1/2 

(3) 



1 + 



2EL 2 



where E and L are the total orbital energy and angular mo- 
mentum, fi is the reduced mass, G the gravitational constant 
and rrih and m 3 are the masses of the halo and the satellite, 
respectively. The second parameter, the pericenter distance 
corresponds to the fictitious closest approach of the objects 
(pretending the system would continue to behave like a per- 
fect two body system long enough) and is given by 




Figure 1. Top panel: cumulative redshift distribution of major 
mergers. The vertical axis shows the fraction of mergers later than 
redshift z in our simulations. Lower panel: Merger rate at redshift 
2, i.e. the number of mergers per unit time normalized to the total 
number of halos with masses above 5.35 X lO n M . 
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We find that 90% of all mergers have bound orbits 
(e < 1) and 71% have values of 0.9 < e < 1. This clear 
maximum of p(e) close to unity can be considered as the 
main feature of the eccentricity distribution. Note that ini- 
tially unbound hyperbolic orbits are as well compatible with 
merger events since orbital energy can be transformed to in- 
ternal one during the merger. 

For the pericenter distance -R pcr i we find a probability 
distribution which can be well described by P(R pcl i/ R v i r ) = 
^exp(-a x i? P cri/-Rvir), where (A, a) « (0.089,5.0) for 
Rpcri/Rvir > 0.2 and (0.12,6.3) for R pcli /R vir < 0.2. The 
mean relative pericenter distance turns out to be 0.18i? v ir- 

We finally checked the distribution of the angles j3 be- 
tween the spin planes of halos and infalling satellites. In 
agreement with expectations we find a distribution that is 
well fitted by dN/d/3 oc sin/3. 

Khochfar & Burkert (2003, KB2003) investigated in 
very detail the statistics of various merger orbital param- 
eters using a huge data sample provided by the VIRGO 
consortium. Apart from minor deviations we find the results 
from this section to be in good agreement with KB2003. 



2004 for a reviews). As will be discussed in this paper, the 
evolution of a dark halo's spin parameter seems as well to be 
a matter of the mass growth history (see also D'Onghia & 
Burkert 2004). Therefore it is worth to illustrate the mean 
merger history of our sample of dark halos throughout their 
evolution. The top panel of Fig^shows the probability for a 
dark halo to have a (final) MM at redshift smaller than z (if 
it has a MM at all), while the lower one shows the number 
of (final) MMs per time, normalized to the total number 
of halos with masses higher than m m i n > 5.35 x 1O 11 M0 
(corresponding to 500 particles minimum). The so defined 
merging rate shows a broad maximum between redshift 1.2 
and 1.8, where the rate is roughly constant. After z w 1 as 
well as before z m 2 the rate is continuously decreasing to a 
value of « 0.02 ± O.OOA/Nh/Gyr and 0.05 ± 0.005/N h /Gyr 
at z = and z — 2.8, respectively, where Nh is the total 
number of halos above m m i n at the corresponding redshift. 
Roughly one third of all MMs occur in redshift intervalls of 
z = — 0.5, 0.5 — 1 and 1 — oo (see top panel of Fig0. This 
shows that the spin parameter of dark halos may undergo 
changes throughout it's whole lifetime, if the spin is affected 
by these merger events. 



3.2 Merger Statistics 

The internal structure of dark matter halos and their embed- 
ded galaxies is to a high degree sensitive to their merging 
histories. Strong evidence exists now that elliptical galax- 
ies originate from MM events (see Toomre & Toomre 1972, 
Naab, Burkert & Hernquist 1999, see Burkert & Naab 2003, 



4 SPIN TRANSFER IN ACCRETION AND 
MERGER EVENTS 

In the commonly established scenario, cosmic structures ac- 
quire angular momentum by tidal interaction with the sur- 
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rounding matter distribution (Barnes & Efstathiou 1987, 
Padmanabhan 1993 and first Hoyle 1949). This mechanism 
works most effectively during the early expansion and col- 
lapse phase, when the large spatial extension of the pro- 
tohalo offers a large lever for tidal forces. The Zel'dovich 
approximation is a powerful tool for an analytical treat- 
ment of the spin acquirement and it turns out that angular 
momentum is growing linearly with time until the point of 
turnaround is reached. 

More recently, cosmologists became aware of merger 
events to have a significant impact on the halos' spin by the 
transfer of orbital to internal angular momentum (Gardner 
2001, Vitvitska et al. 2001 (VK01)). It appears plausible 
that the quantitative evolution of spin transfer in mergers 
is depending on the corresponding orbital parameters. This 
complicates the investigation of the spin evolution in so far 
as we have to deal with nonlinear physics which rules out an 
ab initio analytical description like the linear tidal torque 
theory. 



4.1 Spin Transfer Statistics 

To explore the merger impact on the halo spin parameter in 
detail, the most basic investigation is to analyse the proba- 
bility of a single merger event to rise (or decrease) the spin 
of a halo by a certain amount. The alteration of spin is quan- 
tified simply by the spin transfer 



T\ '.— Afln/Ainit, 



(5) 



where Afi n is the (final) spin parameter of the merger rem- 
nant and the initial Ai n i t corresponds to its main progenitor. 
We want to highlight the characteristics and differences of 
both established versions of dimensionless spin parameter 
definitions throughout this paper. The earlier classical ver- 
sion introduced by Peebles (1980) reads 



A = 



jVe 



/2" 



Here J is the modulus of the total angular momentum, E 
is the absolute value of the total (kinetic plus potential) 
energy and M is the mass of the halo under consideration. 
All quantities are calculated using the set of particles within 
the virial radius J? v j r . Bullock et al. (2001) later defined an 
alternative version, namely 



A' = 



JV4\rJj -^^vir ^vir 



(7) 



where V v i r = GM V1I R~^ is the circular velocity at the 
virial radius R v i r . This latter (dashed) quantity is more fre- 
quently applied in practical analysis due to its independence 
of the density profile. Both quantities are nevertheless ap- 
proximately identical for halos with NFW profiles (Navarro, 
Frenk & White 1996) and concentration parameter c ~ 10 
as well as isothermal halos truncated at R v i r (Mailer, Dekel 
& Somerville 2002). 

Fig. 2 shows the A-distribution for all halos at z — and 
for the merger and accretion dominated subsamples of halos 
as well. One can clearly see a shift between each sample cor- 
responding to the according spin transfer behaviour: As to 
expect, MM dominated halos are shifted to higher A values, 
mM halos occupy lower values. In each case the distribution 
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Figure 2. Distribution of the (classical) spin parameter A. Each 
sample (halos with MMs, without MMs and full sample) is fol- 
lowing roughly a lognormal distribution (see e.g. Cole & Lacey, 
96). Depending on the merger history the distribution is shifted 
to higher (MM sample) or lower (mM sample) values. 
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o-oT \, MM 



(6) 4 




Figure 3. Distribution of the spin transfer rates T\ = Ag n /Ai n i t 
(full lines) and T' x = A~ „/A^j t (dashed lines) for major (circles) 
and minor mergers (squares). Clearly major mergers provide a 
source of angular momentum: roughly two third of all MMs in- 
crease the spin parameter; in the case of minor mergers less than 
half of all events increase A by a moderate amount. Note that 
the behaviours of T\ and are hardly distinct for either merger 
type. See the discussion in the text for details. 



is roughly lognormal. Here we do not show the distribution 
of A' but note that we find a similar behaviour. 

FigE] shows the distribution of the spin transfer rates 
T\ and T' x := Afl n /A[ nit for the case of major (m s > 0.2m^) 
and minor merger/accretion (ra s < 0.2mi,) events. In the 
following discussion all numbers concerning A' are listed in 
brackets behind the numbers corresponding to the undashed 
A. 
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Fig|3| clearly shows that MMs rise the spin parameter 
by a significant amount: 68% (67%) of all MMs rise A (/) . 
In more than 50% of all MMs, A''' is increased by a factor 
of 1.30 (1.31). This is in good agreement with the previous 
study of VK01 who found a median of 1.25 for the dashed 
quantity T' x . For the according mean values we find T\ are 
1.32 (1.34). 

As expected the situation is different for the case of mM 
and accretion events. Again there are both, spin-rising and 
decreasing events, but now with a clear shift to a lower T\ 
regime: The distribution of T\ in this case is significantly 
narrower and approximately symmetric around unity, with 
a linear median value of T\ — 0.98(0.99) and a mean of 1.007 
(0.997). 53% (56%) of all events decrease the spin parameter. 

It turnes out that the behaviour of both spin parame- 
ters is very similar in the above context, regardless of the 
density profile dependence of A/ A' which should lead to vis- 
ible differences especially during the merging phase, where 
the density distribution is strongly time dependent. Never- 
theless the tiny difference of the T\ and T\ distribution in 
the accretion case has striking effects on the global time evo- 
lution of A and A'. We will discuss this important issue in 
detail in S0 
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Figure 4. Cumulative (full line) and differential (dashed line) 
distribution of the ratios of orbital to internal angular momentum. 
In 50% of all mergers L or i, it > 6 X J pr g ■ In less than 4% of all cases 



The diagram demonstrates that the spin transfer 



is dominated by the orbital angular momentum. 



4.2 Dependence on Orbital Parameters 

As mentioned before, the plausible explanation for the quali- 
tative behavior just described is that orbital angular momen- 
tum and the spin of the satellite are transferred to internal 
spin of the remnant during the merging process, according 
to 



T x oc 



J rc 



(8) 



where the indices represent remnant, progenitor, satellite 
and the system's orbit respectively. Exact equality is un- 
likely since some dark matter particles are lost during the 
merging process. is more or less sensitive to the re- 
lative orientation of the vectors on the right hand side of 
Equ. JSJ , depending on the relation of their absolute values. 
Since we are discussing major merger events it is justified 
to assume satellites and halos to have internal angular mo- 
menta of the same order of magnitude. From Fig^Jone can 
see that only in a small fraction (4%) of all MMs the halo 
internal angular momentum is larger than L or bit- In 80% of 
all events -Lorbit > 3 x J prg and still for every third case 
iorbit > 7.7 x J P rg. This demonstrates directly that the spin 
transfer is dominated by the orbital angular momentum. We 
find a small but detectable dependence of T\ on the vector 



angles a^j and a j 



, The bottom panel of FiglBlshows 



a linear fit through the median T\— angle plot. Though suf- 
fering from a large scatter, there is a clearly visible trend 
of increasing spin transfer with decreasing spin-orbit- and 
spin-spin-angle. 

If the transfer of angular momentum is dominated by 
iorbit we should also expect a more elementary dependence 
of T\ on the pericenter distance since (for fixed eccentricity 
e) iipcri is proportional to the satellites' impact parameter 
Pi which, for its part, satisfies L or tu oc vpi. In Fig[51 (upper 
panel) T\ is plotted against R per i/R200 which is well fitted 
by the power law 



R D 



R20 



(9) 



over a R per i/ R200 range of almost 3 orders of magnitude. 
This clearly reflects the leading role of the impact param- 
eter for the amount of spin transfer in a merger event. We 
also investigated the according behaviour of the dashed spin 
transfer T' x and didn't find any significant difference. 



5 BEFORE, DURING AND AFTER MERGING 
— RISING AND DECLINING SPIN 

Major mergers at first appear to be the main sources for 
rising spin parameters in the nonlinear evolution of dark 
halos. The upper panel of FiglrH however offers a more com- 
plex view: it shows the evolution of the mean A along the 
"fictitous" time axis t* , where the latter variable represents 
an individual time scale for each halo with the zero point 
t* — shifted to the halo's final major merger event. In 
other words, the individual curve of each halo j is 

shifted along the (real) time axis such that its final MM 
happens at the same fictitious time t* = 0. The average is 
taken within time intervalls of 0.6 Gyrs. 

For t* < the mean spin parameter is almost con- 
stant or at most slightly increasing with a rate of dX/dt ~ 
(4.8 ±0.8) x 10 _4 Gyr _1 . Within the first 0.6 Gyrs after the 
merger onset A rises fast by 30%, followed by a 5 Gyrs last- 
ing 1st phase during which the spin decreases again with a 
power law A(t*) oc , and a 2nd phase with a further 

decline, following a quadratic function. 

The A'(i*) curve looks similar. There is a significant 
peak after the merger onset and a modest decline later on. 
Note however that both spin parameters have equal mean 
values only at late fictitous times t* , i.e. well after the final 
mergers. For times t* < 3 Gyr A is significantly lower than 
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Figure 5. Top: Spin transfer T\ vs. relative pericenter dis- 
tance. The diagram shows a tight correlation following T\ oc 
(^pcri/^vir) 018 - The error bars reflect the Gaussian tolerance 
for the mean values rather than the standard deviation of the in- 
dividual merger data. Bottom: T\ vs. spin-spin and spin-orbital 
angular momentum angle (ot£j and ajf respectively). As typi- 
cally L or t>it/ Jprg > 1 the T\— angle correlation is clearly weaker 
than the relation between T\ and iJ nr ,-; 




Figure 6. Top panel: Evolution of the mean spin parameters A 
and A . The time axis t* denotes the time for each individual halo 
relative to its final major merger event at t* = 0. The sudden rise 
of A^ ' right after the merger onset is followed by a less rapid 
decline. The middle and bottom panel show the (relative) time 
evolution of the quantities that make up the spin parameters (see 
Eqs. JHJ and Q). In both cases, the sudden rise of the angular 
momentum (open circles) is responsible for the fast increase of 
the spin parameter after t* = 0. Note also the dotted-dashed line 
in the bottom panel, which illustrates the differential growth of 
the virial radius. See text for further discussion. 



A', and the difference increases towards earlier times. For 
t* < A' is also roughly constant. 

To understand the behaviour of the mean spin param- 
eters a' ' in further detail, it is instructive to explore the 
time evolution of the quantities which constitute the spin 
parameters, namely the angular momentum L(t*), the to- 
tal energy E(t*) and the total mass M(t*) for the case of 
A (Equ. JSJ) and the virial radius R V [ T and velocity V v i T for 
A' (Equ.0. The results are shown in the middle and lower 
panel of Fig|S| Each curve is normalized to unity at the 
merger onset t* = 0. 

We consider the middle panel of Fig|S] first, represent- 
ing the undashed mean A. Already for t* < the perma- 
nent smooth matter infall causes a steady growth of the 
total energy absolute value \E\, angular momentum L and 
mass M within the redshift dependent virial radius. Until 
the merger's onset the spin rising effect of \[E and L com- 
pensates and somewhat exceeds the decreasing effect of the 
M -5 / 2 term in the denominator of Equ|S| While the values 
of L and \J~E rise by factors of approximately 5 (neglecting 
the initial rise which is likely due to the poor statistics for 
t* < —8 Gyrs) and 2 respectively, M -5 ^ 2 drops by a factor 
of m 7 during the same period. 

At t* — 0, a large satellite enters the halo and the mean 



angular momentum and the total energy rise by factors of 
2.3 and 1.4 respectively within only 0.6 Gyrs, while the total 
mass term (M -5 / 2 ) is decreasing by a factor of « 0.25. One 
must be carefull not to simply multiply these quantities in 
order to predict the beh aviour of A. N ote that we are con- 
sidering mean values (L, M~ 5 / 2 , E 1 / 2 ) here, i.e. the sums of 
single values, while the mean spin parameter X(t*) in Figlulis 
a sum of the products of individual y^Ej, Lj and M~ 2 ' 5 val- 
ues. Both quantities, A(t*) and L(t) x E 1 / 2 ^*) x M~ 2 - 5 (t*) 
would be comparable only for small variances of the single 
values. 

We can nevertheless draw some qualitative conclusions 
on the roles of L, E and M for the evolution of the spin 
parameter X(t*). Obviously, the rise of angular momentum 
is the quantity mainly responsible for the sudden rise of A 
after t* = while the total energy seems to contribute to 
a minor degree. After t* ~ 0.6 Gyr, when A is decreasing, 
the figure shows that this is due to the change of the mean 
mass, which is continuing to increase by a factor of ~ 5 due 
to accretion between the final merger and z = while L and 
E 1/>2 are not changing significantly during the same period. 

The lower panel of Fig.© shows the situation for the 
case of A': Before t* — the increasing mean mass and virial 
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6.1 Global Time Evolution 




Figure 7. Distribution of the virial coefficient r\ at various red- 
shifts. Note that even at z = 0, the distribution doesn't peak 
at r] = 1, but is shifted to slightly higher values due to infalling 
material. 



radius is almost canceling the effect of the increasing angu- 
lar momentum. At t* — the spin parameter is rising by 
more than 30% within 0.6 Gyrs, caused by the fast rising 
angular momentum as in the (undashed) A case. Later on, 
i.e. after t* = 0.6 Gyr the growing mass and radius in com- 
bination with an almost constant mean virial velocity and 
angular momentum of the halo sample is responsible for A 
to decrease. 



6 THE VIRIAL COEFFICIENT rj 

We store our simulation data at 21 particular times between 
z ~ 2.5 and z = with approximately equal time intervalls 
At ~ 6.3 x 10 8 yrs. We therefore can detect a MM in a phase 
within a period At after the individual virial radii of both 
merger partners had first contact. In such an early stage 
of merging the merger remnant is very disturbed and far 
from dynamical equilibrium. This might cause the tempo- 
rary peak of the spin parameter values at t* ~ 0. To inves- 
tigate this question and quantify the relaxation of a halo we 
employ the virial parameter 



77 := 



2-Ekin 



ipot| 



(10) 



where i?kin and E po t are the total kinetic and potential ener- 
gies of all particles within R v i T . According to the virial theo- 
rem rj should approach unity for an isolated relaxed object. 
Note that infalling material is affecting the surface pressure 
S p = — J pr ■ dS of dark halos in cosmological simulations (p 
is the pressure with radial direction f on the surface element 
dS). Since the scalar virial theorem in absence of magnetic 
fields reads 2E^ n + E pot + S p = 0, a momentum flux induced 
negative surface pressure causes 2Eki n /\E po t\ > 1 (Shapiro 
et al. 2004), which we indeed find for the mean virial co- 
efficient rj at any redshift (Fig. 8). Figure 7 shows the dis- 
tribution of rj and its dependence on redshift. Only a small 
fraction of halos have r\ = 1. In addition, rj on average in- 
creases with increasing z while the width of the distribution 
remains roughly constant. 



Before we are going to consider the correlation between the 
spin parameters and the virial coefficient, we want to ex- 
amine the time and redshift evolution of rj. As shown in 
the top panel of Fig|H|the mean virial coefficient is clearly 
decreasing with decreasing redshift, well described by the 



quadratic function r](z) 



-3.3 x 10-\z + 2.7Y + 1.3; this 



corresponds in good approximation to a linear decline of rj 
with time, namely rj(t) w 1.35 — 1.7 x 10~ 2 i/Gyr. It reflects 
the general tendency of halos to approach a more dynami- 
cally relaxed state, corresponding to rj fa 1, according to the 
fact that interactions between halos in later times are less 
frequent and disturb halos less than at high 2. 
In addition, the surface term S p becomes less important. 
The lower panel of Fig|5] shows the evolution of rj with fic- 
titous time t* relative to the final merger (analogous to the 
spin parameter evolution in O . We find again the peak- 
like feature at t* = 0, similar to X^''(t*). At the epoch of 
the major merger happens, the general decline stops and rj 
increases temporarily instead. Compared to the merger in- 
duced rise of A(t*), the change in rj is measurable already 
before the merger onset: More than 2 Gyrs before t* = 0, 
rj(t*) starts to deviate significantly from the linear decline 
which is shown by the dashed curve in Fig. 8 and already 
0.6 Gyrs before merging rj experiences its fastest rise. Ob- 
viously the internal dynamical equilibrium of the halo is 
already disturbed by the proximity of a satellite; other than 
in the spin parameter case, the effect on rj is preceding the 
physical contact, rj reaches it's maximum value at the same 
time as A. About 2.5 Gyrs after the merger, rj is continuing 
to evolve as if nothing happened, i.e. rj(t*) is following the 
same linear behaviour as before (see Fig|5Jl. Obviously the 
temporary fast increase in the spin parameters shortly after 
a major merger is at least partly caused by the fact that the 
halos are out of virial equilibrium and the spin parameters 
achieve again a more typical value after a relaxation time 
which corresponds to a few dynamical timescales. 



6.2 Correlation between rj and 

Peebles' classical definition of the spin parameter (Equ|SJ 
makes use of both the kinetic and potential energy. There- 
fore one can expect a direct correlation between A and 
rj — 2-Eidn/ |-Epot I, disturbed only by the scatter of the inde- 
pendent angular momentum J in the definition of the spin 
parameter. The Bullock et al. alternative A' is not using me- 
chanical energies explicitly. It should nevertheless depend on 
relaxation since the virial radius which defines the region in- 
side which A is determined, is not well defined in a disturbed, 
unrelaxed halo. 

Fig[§] shows the correlation between A^'' and the virial 
coefficient rj. Both A and A' are increasing for higher values 
of rj, following 



A 



(') 



a + (3rj 



within a range of 0.9 < rj < 1.5, which contains more than 
90% of all halos. Two clear trends are visible: 

(i) The Bullock spin parameter A' is significantly stronger 
correlated with the virial coefficient rj. 
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Figure 8. Top panel: Redshift evolution of the mean virial coef- 
ficient, rj is decreasing with declining redshift, reflecting the ten- 
dency of halos to become more relaxed and infall to be less impor- 
tant. The quadratic decline with decreasing redshift corresponds 
approximately to a linear decline with time. Lower panel: The 
linear behaviour of rj with victitious time t* is disturbed only by 
the merger event, analogous to the behaviour of the spin param- 
eter (grey line). Before and after the merger rj(t") is described by 
the same linear function (dashed line, see text for details). 



(ii) The correlation becomes slightly weaker at higher red- 
shifts. 

The values of the parameters a and /3 can be found in 
the caption of Fig[5] 

We finally examine the ratio of the spin parameters A'/A 
in order to further illuminate the differences in their be- 
haviours. Fig llUI shows that A'/A is decreasing with time - 
in quite a similar manner as the virial coefficient - again 
featuring a hump in the vicinity of t* = 0. The strong corre- 
lation between A'/A and rj is illustrated by the rj/{\' /A) plot 
in the same figure: the linear behaviour of this quantity is 
only disturbed for a relatively short period directly after the 
merger onset, indicating that A'/A is falling back to a lin- 
ear evolution more quickly than the virial coefficient. Note 
that r]/(\'/\) in Fig. HUl is not yet affected at t* — when 
the individual quantities already reach values close to their 
maximum ones. Thus the ratio A'/A is stronger correlated 
with the virial coefficient rj than each spin parameter sep- 
arately. This result turns out to be important for the next 
chapter. 




4 



Figure 9. Spin parameter dependence of the virial coefficient rj. 
Each case can be well fitted by A^ = a + (3 X 10 -3 ?? 4 . The fit 
parameter values are (a, 0) = (0.036, 2.4) for A; z > 0, (0.028, 7.9) 
for A'; 2 > 0, (0.037, 4.7) for A; z = and (0.028, 11) for A'; z = 0. 
A''' {v = 0) is almost independent of the redshift. Note that "z > 
0" refers to the complete sample of halos excluding only those at 
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Figure 10. Evolution of the spin parameter ratio A'/A (full 
squares) with respect to the final merger event (t* = 0). Plot- 
ted are also the virial coefficient rj (open circles) and the ratio 
?7/(A'/A) to illustrate the strong correlation between r) and A'/A, 
which is only disturbed for a short time period after the major 
merger. 

7 TIME EVOLUTION OF THE SPIN 
PARAMETERS 

In this chapter we want to address the question, how the 
sum of merging and accretion affects the global evolution of 
the spin parameters with cosmological time. It turned out in 
this work and previous ones (e.g. Vitvitska et al. 2001), that 



individual MM and mM events can either rise or decrease 
the spin parameter: While mMs give T\ ~ on average, it 
is obvious that MMs tend to clearly rise A. 

The spin history of any single halo is highly chaotic as 
soon as the smooth phase of spin- up increase by tidal torques 
is left behind and non linear matter acquisition comes into 
play (Peirani, Mohayaee & de Freitas Pacheco 2003, Primack 
2004) . But if we instead consider the mean value of all halos 
we may still expect a smooth redshift evolution of the mean 
spin parameter due to the cumulative increase of the total 
number of mergers during the cosmic evolution shown in 
FigEJ 

We find this to be valid for both spin parameters with a 
nevertheless important and surprising difference. In Fig llll 
we see that A , averaged over all halos is approximately con- 
stant between z — and z ~ 2. The mean value of the 
undashed classical A in contrast is significantly increasing 
with cosmological time, following approximately the linear 
relation 



A(z) = 0.039 - 3.5 x 10" 3 2. 



(11) 



It is instructive to consider just the mM/accretion dom- 
inated subsample of halos. First we note that the values of 
both A and A' lie well below the full sample of (MM+mM) 
mean values as a result of the missing angular momentum 
"kicks" by major mergers. But while A due to mMs de- 
creases significantly with decreasing redshift, we find that 
the opposite is true for the undashed classical A. This un- 
veils the crucial difference between both spin parameters. 

To understand this feature we consider again FigEl 
which shows the spin transfer distributions. As already men- 
tioned, one would hardly speak of distinct behaviours be- 
tween the transfer rates at first glance: We noted values of 
T\ = 1.007 and T' x = 0.997 for the mM halos. In our sim- 
ulation every halo suffers this rise (or reduction) of its spin 
parameter around 18-19 times 2 on average between z — 2 
and 2 = 0. The net effect of all individual mM events is thus 
AA (,) := A^ /AW a = Pf'] 19 • For the dashed version this 
gives AA' = 0.997 19 « 0.94, while for the classic version we 
get AA = 1.007 19 « 1.14. FigHTl shows that these numbers 
correspond roughly to the redshift evolution of either minor 
merger sample! Thus the tiny differences between the T\ and 
T' x distributions lead to a distinct time evolution of the spin 
parameters. It is beyond the scope of this paper to analyse 
the detailed physical mechanisms for this difference. 

According to the Fig llll major mergers in either case 
have only corrective character: Though individual MMs rise 
the spin parameter by 130% on average, they occur too 
rarely to overtrump the net effect of accretion and minor 
mergers which act quasi permanently. We conclude that mi- 
nor merger events and accretion rather than major mergers 
provide the driving force guiding the behaviour of the spin 
parameters. Their very distinct time evolution is entirely 
caused by tiny differences of the spin parameters' response 
to accretion, while major mergers yield a minor additional 



2 Remember that we defined minor mergers/accretion events as 
any mass acquisition between two data outputs that can not be 
classified as major merger. The latter happens about once during 
a typical halo life. Since we have 20 output dumps since z = 2, in 
19 of them we have to deal with mMs by definition 
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Figure 11. Redshift evolution of the mean spin parameters A and 
A for the mM- and the full sample of halos. While A is increasing 
with time for both halo samples, we find A const for the full 
halo sample and A decreasing for the accretion dominated halos 
in agreement with VK01. 



effect, acting similarly on either A and A'. Minor mergers 
continuously bring in a net amount of angular momentum. 
Large scale gravitational torques therefore seem to dominate 
the origin of angular momentum in galaxies in conflict with 
previous ideas that continuous infall occurs from random 
directions, leading to a zero net effect on A™. 



8 SUMMARY AND DISCUSSION 

We performed simulations of large scale structure formation 
in a ACDM universe and extracted all Friends-of-Friends 
selected dark matter halos at 20 distinct redshifts between 
z = 2.8 and z = 0. Our principal goal was to investigate the 
differences of the two commonly used versions of the spin 
parameter, to understand their behaviour under the influ- 
ence of halo major mergers and accretion and to determine 
whether the spin parameter evolution is driven by minor ao 
major mergers. 

The spin transfer T^' — Ag^/A^L for any minor (mM) 
and major merger (MM) event throughout the simulations 
was calculated . We find that (1) the probability distribu- 
tion T'x agrees well with the one found by Vitvitska et al. 
(2004) and (2) the corresponding "classical" quantity T\ is 
distributed almost identically. This holds for either minor 
and major mergers regardless of the considerable distinct 
values of A and A' especially during merger events. 

We analysed the spin transfer dependence of the orbital 
parameters. The pericenter distance between a halo and an 
infalling satellite was identified as the free parameter which 
is mainly responsible for the spin transfer . This is a re- 
sult of the high values of L orb /Ji nt > 3(5.7, 11) in 80% (50%, 
20%) of all major mergers (L or b and Ji n t being the orbital 
and internal angular momentum, respectively). It is mainly 
orbital angular momentum that is transferred to internal 
spin during mergers. The relative orientations between halo 
and satellite spin vectors turned out to have a minor (but 
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detectable) influence on , resulting from similar order-of- 
magnitude values of internal and orbital angular momentum 
in a considerable amount of merger events. 

In order to further analyse the effect of merger events 
on the evolution of the spin parameters, we plotted the evo- 
lution of A''' along a fictitious time axis which describes 
the time relative to the final MM event of each individual 
halo. It turned out that after the merger onset, A^ expe- 
riences a peaked rise of ~ 30%, and a considerable decline 
later on. The first is due to relaxation effects, combined with 
the transformation of orbital to internal angular momentum, 
while the latter seems to be the result of subsequent mass 
accretion in combination with a constant mean value of the 
added specific angular momentum. 

The time evolution of the virial coefficient rj = 
2-E]dn/|-Epot| and its correlation with the spin parameters 
was investigated. A''' is very sensitive to rj following A 1 -'' « 
a + /3rj 4 . The values of a and /3 depend on redshift and 
the type of the spin parameter. The correlation is stronger 
for lower redshifts and for the dashed spin parameter A'. 
On the fictitious time axis t*, rj shows the same peak-like 
feature around t* = as the spin parameters, though less 
pronounced. A strong correlation was found between the ra- 
tio A'/A and rj. While one might not be surprised about a 
A (,) — rj correlation in principal (due to the direct or indirect 
use of the mechanical energies) the detailed origin of this 
particular behaviour is still an open issue. 

We finally examined the redshift evolution of the 
spin parameters. While the global mean value A turnes 
out to be roughly redshift independent, we find that A is 
significantly increasing with time. If we restrict the halo 
sample to those which never had major mergers, the picture 
changes slightly: While A' is still significantly growing with 
time, we find A to decrease moderately. The earlier result of 
similar distributions is not contradicting this result: In 
fact A and A' have mean spin transfers which are different 
by a small but important amount for the mM sample 
of halos (T^ = 1.007(0.997)). This causes a significant 
difference in the global evolution since the "operation" 
A^ — > x \W is acting about 18 times on each halo 
between z = 2.8 and z — 0. The distinct behaviour of 
both spin parameters is thus induced by their distinct 
reactions on minor merger and accretion events. Due to 
their comparable rare appearance major mergers must be 
considered as only minor additional corrections to the A 
evolution which is dominated by a coordinated action of 
minor mergers caused probably by large scale stream flows 
that could be the result of tidal torques. 

We like to thank Hans- Walter Rix for enriching com- 
ments and general support. This work was supported by the 
Deutsche Forschungsgemeinschaft via SFB 439, Galaxien im 
jungen Universum. 
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